Dynamics of domain growth driven by dipolar interactions in a perpendicularly 

magnetized ultrathin film 
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Measurements of the ac magnetic susceptibility of perpendicularly magnetized Fe/2 ML 
Ni/W(110) ultrathin films show a clear signature of the dynamics of domain growth and domain 
density changes in the striped domain pattern that this system supports. The susceptibility peak 
measured at different constant heating rates in the range 0.20 K/s < R < 0.70 K/s shifts to higher 
temperature as the heating rate is increased. Analysis using a relaxation model demonstrates quan- 
titatively that the dynamics is driven by a non-equilibrium domain density at (nearly) zero field 
(i.e. by dipole interactions), and that the temperature shift is due to a response time determined 
by the pinning of local domain wall segments by structural defects. The fundamental time scale for 
relaxation of the domain density driven by dipole interactions is of order 10 times slower than the 
fundamental time scale for an individual Barkhausen step driven by an applied field. The increase 
in the fundamental time scale reflects the relative size of dipole and Zeeman energies, and the need 
for the correlated motion of the many local domain wall segments required to affect domain growth. 

PACS numbers: 



I. INTRODUCTION 

Ultrathin films with perpendicular anisotropy are an 
example of a larger class of two dimensional systems 
where strong, short range attractive interactions and 
weak, long range dipole repulsive interactions lead to the 
formation of domain patterns*^— The magnetic domain 
patterns in these films provide a unique opportunity for 
the study of one-dimensional domain walls (more prop- 
erly domain "lines" ) in a two dimensional magnetic sys- 
tem and to understand the way in which domain wall 
dynamics and fluctuations determine many of the finite 
temperature magnetic properties of the film. 

The present article outlines a quantitative experimen- 
tal investigation of the dynamics of magnetic domain 
growth driven exclusively by long-range dipole interac- 
tions. It is, perhaps, surprising that, given the intense 
interest in these systems over the last few decades, this 
question has not been addressed quantitatively. The ab- 
sence of clean experiments is due to the complicated 
hierarchy of magnetic relaxation processes, spanning a 
wide range of time scales, that control the dynamics of 
these systems. At sufficiently low temperature, the films 
support an ordered stripe domain pattern^ The motion 
of these existing domain walls, when driven by an ap- 
plied magnetic field, is controlled by localized pinning of 
the domain walls at microstructural defects in the film. 5 
As the temperature of the films is increased, the stripe 
domain density changes exponentially with temperature 
through domain growth and creation mechanisms that 
are driven by long range dipole interactions Finally, 
at a sufficiently high temperature domain wall fluctua- 
tions can drive a transition to a different domain pat- 
tern through the proliferation of topological defects in 
the pattern itself.— In order to study the domain growth 
dynamics driven by dipole interactions, it is necessary 
to use a zero (or very low) field technique with a very 



wide dynamic range. Only then can one identify the en- 
tire heirarchy of dynamic processes, isolate the desired 
mechanism and study it on the appropriate time scale. 

Almost all of the existing experimental studies do 
not meet these requirements because they have used ei- 
ther large magnetic fields or static imaging. Dynami- 
cal studies have used large applied magnetic fields ei- 
ther to study a magnetically saturated state very far 
from equilibrium, 10 or to move domain walls using Zee- 
man magnetic forces rather than weak internal dipole 
forces i 11 ' 12 Studies of the domain structure often use 
magnetic microscopy in no applied field, but they are lim- 
ited to a near-static characterization of the domainsi^i^ A 
very few studies infer dynamical processes qualitatively 
by noting the disappearance of individual domains in suc- 
cessive static domain images^^ but most recognize the 
dynamic limitations of magnetic microspcopy by relying 
on static imaging of domain changes as a function of film 
thickness^— rather than as a function of temperature. 

In contrast, measurements of the ac magnetic suscep- 
tibility cannot directly observe domain geometry, but do 
offer access a wide range of time scales using very small 
magnetic fields that do not overwhelm the dipolar in- 
teractions. There are a handful of reports that use ac 
susceptibility to characterize domain wall dynamics. On 
a short time scale, a small ac field is used to "wiggle" 
existing domain walls in order to study domain wall pin- 
ning by microstructural defects 5 -*^ in the film that act 
on a fundamental time of 10 -9 s. On a long time scale, 
measurements of the susceptibility curve at different slow 
heating rates, R < 0.1 K/s (so that the entire curve is 
traced in many minutes to hours), have revealed the dy- 
namic processes involved in resolving topological defects 
in the domain pattern itself as the system undergoes a 
phase transition from one pattern to another— These 
collective processes act on a fundamental time of 10° s. 
The present article reports experiments on Fe/2.0 ML 
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Ni/W(110) films, and concentrates on the intermediate 
time scale, where dipole interactions drive changes in the 
domain density as the temperature is increased. Mea- 
surements of the ac susceptibility curve as function of 
faster heating rates, R > 0.1 K/s, reveal that the peak 
of the susceptibility shifts to higher temperature with 
greater R. An activated model of domain growth driven 
by the non-equilibrium domain density at zero field re- 
produces these results over a range of heating rates and 
sample thicknesses using a single value of a single ad- 
justable parameter. This parameter, the characteristic 
time scale for dipole-driven domain growth in these films, 
indicates a dynamics which is ~ 10 5 slower than that 
for the field-driven motion of existing domain walls, but 
much faster than that involved in removing pattern de- 
fects. Thus, the use of ac susceptibility as a probe of the 
domain dynamics permits a coherent description of the 
entire hierarchy of competing processes which determine 
the evolution of these systems. 

The remainder of the article is divided into 4 sec- 
tions. In section II, quantitative models that describe 
the domain wall dynamics on different time scales are 
reviewed and/or developed. Section III presents exper- 
imental results of ac susceptibility measurements on ul- 
trathin Fe/2.0 ML Ni/W(110) films, and these results 
are analysed quantitively in section IV. The final section 
summarizes the findings of these studies. 



II. THEORY 



Local movement and pinning of existing domain 
walls: t„ 



is increased, thermal fluctuations renormalize the mean 
field anisotropy and A(T) is reduced. This causes the 
domain wall energy to decrease and the domain density 
increases exponentially with temperature. 

The magnetic susceptibility of a domain phase is due 
to the motion of existing domain walls when a small field 
applied perpendicular to the surface causes the width of 
domains with a parallel magnetization to grow by S , 
while those with an antiparallel magnetization shrink by 
S. If the average domain density is n , then the resulting 
magnetization 



M = M sat n5, 



(2) 



where M sat is the saturation magnetization. In the limit 
of a small applied field, the equilibrium dc susceptibility 



ndn eq (T) 



Aexp(— kT), 



(3) 



where d is the film thickness, and A and k are phe- 
nomenological parameters that have been shown to de- 
scribe the effect of the exponential increase in domain 
density in experimental susceptibility data^^ 

The ac susceptibility measures the oscillation of the 
domain walls in a small sinusoidal field at angular fre- 
quency ui. The response is retarded by pinning of micro- 
scopic sections of the domain walls at structural defects 
with an average binding energy E p i„. The response time 
T p can be modelled by an Arrhenius Law with 



t p (T) = r 0p exp(E pin /kT). 



(4) 



Ultrathin (1-10 atomic layers thick) ferromagnetic 
films may have a magnetocrystalline anisotropy that 
favours perpendicular magnetization. In some cases, it is 
strong enough to overcome demagnetization effects, leav- 
ing a small residual perpendicular anisotropy. — The long 
range magnetic dipole interaction in this geometry is an- 
tiferromagnetic and leads to the creation of domains. 17 
The integration of the dipole interaction over a two di- 
mensional film gives a logarithmic dependence of the en- 
ergy on the domain density. Inverting this relation yields 



the equilibrium domain density^ 



18.19 



(T): 



n eq {T) 



2 E W (T) 

— exp( 

TTt llm 



1). 



(1) 



In this equation £ = Wy/T/X(T) is the domain wall width, 
and E\y = Ay/TX(T) is the domain wall energy per unit 
length. In these definitions, F is the domain wall stiff- 
ness derived from the nearest neighbour magnetic ex- 
change coupling, and A(T) is the effective perpendicular 
anisotropy (including contributions from both the crys- 
talline and demagnetization terms), fl is a constant that 
sets the scale of the magnetic dipole energy, and m is the 
thickness of the film in monolayers. As the temperature 



The ac magnetization is then given by a relaxation equa- 
tion with 



dM(t) 
dt 



-1 



(M(t) - X e q (T)H(t)). 



The steady state solution is 



X(T) 



iu>T p (T) 



1 



) o t-a Xeq{T)- 



(5) 



(6) 



The dynamic prefactor in eq.© describes the effect of 
pinning. The susceptibility decreases at high tempera- 
ture due to the increase in domain density, and at low 
temperature due to pinning. A peak situated roughly 
where ujTp(T P i n ) = 1 divides regions where the domain 
walls are pinned or free. The characteristic pinning tem- 
perature is therefore 



T 



-E n 



/eln(wrop)' 



(7) 



Studies have shown that ro p ~ 10 9 s for this microscopic 
pinning process. 
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B. Removal of topological pattern defects: Td 



C. Changes in domain density: r n 



Theoretical and numerica l 20 ' 21 studies have indi- 
cated that a phase transition between the striped domain 
phase and one of a number of delocalized domain phases 
is expected at a temperature sufhcently high that the do- 
main density is large, and domain wall fluctuations are of 
over-riding importance. The delocalized phases are char- 
acterized by the proliferation of topological defects in the 
stripe domain geometry that break the stripes into seg- 
ments, and may re-orient the segments, so that positional 
and/or orientational long range order are lost. 

Numerical simulation s 22 ' 23 have suggested that when 
the system is quenched from a high temperature delo- 
calized phase to the low temperature striped phase, the 
topological defects persist for a very long time. This is ei- 
ther because macroscopic rearrangements are required, 22 
or because the transition is non-continuous ^ A recent 
experimental study of 1.5 ML Fe/2ML Ni/W(110) films 
has observed the relaxation of the topological defects 
in quenched films indirectly, through measurements of 
the ac susceptibility.-— The films were cooled from high 
temperature (360 K) at a rate of i?=-0.10 K/s, and the 
susceptibility curve was measured for different constant 
rates of heating. The whole susceptibility curve was seen 
to shift to higher temperature when the constant heating 
rate R at which it was measured, was decreased. The 
long time scale (many minutes to an hour) on which 
this occurred, Td(T), was identified as the characteris- 
tic time for topological defects to relax. For large R, the 
measurement time was much less than Td(T), and sys- 
tem retained the defects of the delocalized phase, giving 
a susceptibility curve with an instrinsically lower peak 
temperature. For small R, the measurement time was 
many times Td(T) and the system relaxed to the equilib- 
rium, striped phase that had an intrinsic peak at higher 
temperature. Thus the peak temperature decreased with 
increasing R. 

A quantitative description of the relaxation of topolog- 
ical defects was provided by 



T d {T) = T 0d e-xp(E d /kT), 



(8) 



where the fundamental time scale rod was found to be of 
order 10° s, and Ed is the barrier to the removal of the 
topological defect. The peak temperature as a function 
of the heating rate was well described by 



T pk (R) = T - Aexp(-t eff (R)), 



(9) 



with t e f f (R) the effective number of time constants that 
have passed during the measurement while heating at 
rate R from initial temperature Tf. 



teff(R) 



T t Rrd{T) ' 



(10) 



To is the peak temperature when the relaxation to the 
equilibrium stripe phase is complete. 



Since the entire susceptibility curve relaxes along the 
temperature axis with time constant Td, the fundamen- 
tal time scale for changes in domain density driven by 
dipole interactions must lie between the two extremes 
set by tq p and rod- It is possible to access this time scale 
as well, using measurments of the ac susceptibility with 
small applied fields. 

Measuring x{T) involves changing the temperature at 
a rate R (K/s). Because structural defects pin the do- 
main walls, it takes time for domains to grow or contract, 
and the domain wall density n(T) will lag behind the 
equilibrium value n eq (T) by a relaxation time r„(T). In 
this case, the measured susceptibility will be 



,„x 1 — icJTt,(T) , , 



Xeff(T) 



Trdn(T) ' 



(11) 



(12) 



where n(T) is the history-dependent present value of the 
domain density during heating or cooling. It is deter- 
mined by a relaxation equation that is developed in Ap- 
pendix A. The constant rate of temperature change, i?, is 
used to remove the explict time dependence from ea. (|A5[) . 
so that the relaxation equation governing changes in do- 
main density is 



-1 



dn(T) _ 
dT ~ RctT p (T) 



(n(T)-n eq (T)). 



(13) 



An important point in the development of this equa- 
tion is the relation between the time constant for the 
field-induced oscillation of existing domain walls, t p (T), 
and the time constant for dipole-induced changes in the 
domain density, r„(T). Since the pinning sites that re- 
tard the motion of existing domain walls when a magnetic 
field is applied also impede the growth of new stripe do- 
mains driven by dipolar interactions, the activation en- 
ergies for the two processes are the same. However, the 
characteristic times Ton is expected to be different than 
Top because of the need for correlated motion of a larger 
region of domain wall. For instance, whereas the mo- 
tion of any part of an existing domain wall will affect the 
response to a magnetic field equivalently, the motion of 
some regions of the domain wall are much more effective 
than others in changing the domain density. Changing 
the domain density at constant magnetization involves 
a large movement of a small length of domain wall near 
the end of a domain segment. Finally, in the relaxation 
approximation, the number of active growth sites is as- 
sumed to be proportional to the difference between the 
present value of the domain density and the equilibrium 
value, but the value of this proportionality constant is 
not obvious from first principles. All of these geometric 
factors are taken into account through an empirical factor 
a in eq. (|A5[) . so that To n = a.TQ P and t„(T) = otr p (T). 
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It is straighforward to solve eq. (fT3|) using the phe- 
nomenological expansion for Xeq(T) in eq.Q. At high 
temperature the pinning is ineffective, and n(T) — ¥ 
n eq (T), whereas at low temperature the pinning is so 
effective that dn/dT — > and n(T) saturates. T n is a 
characteristic temperature dividing these two regimes. It 
can be estimated^ by setting = 5^ to give 



n(T) = (1 - RaKT p (T))n eq (T). 



(14) 



When RaKT p (T n ) = 1 the model becomes unphysical be- 
cause the domain density cannot change quickly enough 
to maintain equilibrium, and the domain density must 
saturate. This occurs near the temperature 



T, 



-En 



kla(RKaTop) 



The ac susceptibility in eq. (fTTj) depends on the relative 
values of the characteristic temperatures for domain wall 
motion, T p i n , and domain growth, T n . If T n < T p i n , then 
pinning will stop the oscillatory motion of the domain 
walls in the applied ac field, even though the domain 
density can still change through domain growth. The 
dynamical factor involving ujt p will cut off the suscepti- 
bility, so that it is insensitive to changes in the domain 
density below T p i n . In this case, Xeff ~ Xeq over the 
temperature range of the susceptibility peak, and the sus- 
ceptibility will not depend upon R. On the other hand, 
if T n > T p i n i then the saturation of the domain density 
occurs in a temperature range where the field-induced 
oscillation of the domain walls is not pinned, and the 
susceptibility gives a robust signal. The curve shifts to 
higher temperature with increasing R, since increasing R 
increases T n in eq. (fT5)) . (Note that the direction of this 
shift is opposite to that described in the previous section 
for the relaxation of topological pattern defects.) Equat- 
ing T p i n and T n gives the condition that divides these 
two behaviours: 



an 



R' 



(16) 



Choosing measurement parameters on the right side of 
this equation such that they are smaller than the com- 
bination of physical parameters on the left side, ensures 
that peak temperature of the susceptibility, T p k, is sensi- 
tive to the relaxation of the domain density when heating 
at rate R. The experimentally practical range of heat- 
ing rates is from 0.20 K/s (below which the relaxation 
of topological defects dominates^ 5 -), to 1.00 K/s (above 
which the decrease in measurement time at any temper- 
ature results in too much noise). A typical value of k is 
0.04 K~ l . Choosing to — 210 Hz then gives the minimum 
value a m in ~ 10 4 for which the measurement will be sen- 
sitive to the relaxation dynamics of the domain density. 
This magnitude for log 10 (aro p ) is comfortably near the 
middle of the range between the experimentally deter- 
mined values of log 10 (ro p ) and log 10 (Tod). If the actual 
value of a is smaller, then the experimental value of lo 
can be reduced. 



Finally, there is the question of whether or not the 
size of the shift in T p j, as a function of heating rate will 
be large enough to measure. Appendix B presents an 
estimate of the expected size. According to eq. (|B6[) . the 
linear term in the shift of the peak temperature of the 
susceptibility curve, T p k, evaluated locally at R = Rq, is 
given by 



dT 



pk 



dR 



B 



1 



Ro(yp 



3Epin 

kT* k 



+ k) 



(17) 



From the measurements in ref. (fH) for 1.5ML Fe/2 ML 
Ni/W(110) films, eq.(H3 gives B « +5 s for R = 0.5 
K/s. This small slope is consistent with the data, but 
cannot be reliably detected. The parameters E pin , K and 
T p k all decrease quickly with film thickness. 5 Films with 
2ML Fe are estimated to have values of B larger by a 
factor of 2 to 5. This should provide a clear signature of 
the dynamics of domain wall growth. 



III. EXPERIMENTAL RESULTS 

Domain dynamics were studied in a series of Fe/2 ML 
Ni/W(110) films. The nickel layer in this structure estab- 
lishes a slightly strained, epitaxial f.c.c. (Ill) template 
on the tungsten substrate, upon which f.c.c. Fe grows epi- 
taxially for a few layers^ 5 - Due to the f.c.c. Fe/vacuum 
interface, this system exhibits perpendicular anisotropy. 
Above 2.2 ML Fe, the perpendicular anistropy exists only 
at low temperature and the magnetization becomes in- 
plane at a spin reorientation transition^ Below 1.25 ML, 
the magnetic behaviour is complicated by the incomplete 
formation of the Fe layer. In the present study, the film 
thickness was restricted to the intermediate range where 
the description in section II. A has been shown to be valid. 

The films were grown in ultrahigh vacuum using ther- 
mal evaporation from a evaporator with a calibrated in- 
ternal flux monitor, following a procedure established in 
previous studies^ 5 - The first monolayer of Ni was an- 
nealed to 600 K to ensure good wetting and the growth 
was monitored using Auger electron spectroscopy and low 
energy electron diffraction. The substrate temperature 
was measured using a W-Re5%/W-Rel0% thermocouple 
embedded in the W crystal, and the temperature was 
controlled by a combination of static cooling through 
a copper braid attached to a liquid nitrogen reservoir, 
and active heating by radiation from a filament just be- 
neath the crystal. The rate of heating or cooling could 
be kept constant to within 0.05 K/s. The maximum con- 
trolled cooling rate across the whole temperature range 
was -0.10 K/s. The maximum useful heating rate was 
1.0 K/s, above which the reduced number of data points 
introduced excessive noised 

The ac magnetic susceptibility was measured using the 
polar magneto-optic Kerr effect^ where the rotation of 
the plane of polarization of linearly polarized light is pro- 
portional to the perpendicular component of the magne- 



tization. A current coil near the sample created a si- 
nusoidal ac field of amplitude 2.0 Oe, and the very small 
polarization rotation from the ultrathin film was detected 
using a nearly crossed polarizer and lock-in amplification 
of the output from a photodiodc. The ac field had a fre- 
quency of 210 Hz in these studies. The ac current through 
the heating filament had a much higher frequency, so that 
it did not interfere with the measurements. 

Because the dynamics being studied change the shape 
of the susceptibility curve and shift it in temperature, 
considerable effort was invested in developing procedures 
that removed systematic errors and resulted in repro- 
ducible data traces. These included: annealing the films 
to 400K before measurements commenced, and subse- 
quently never heating above 360K; cooling the sample 
from 360 K at a rate of -0.10 K/s before each susceptibil- 
ity measurement was recorded during heating; discarding 
the first heating trace to ensure a consistent magnetic his- 
tory for each measurement; randomizing the order of the 
measurements for different values of heating rate; demon- 
strating that the lock-in time constant of 2 s had no effect 
on the shape or position of the susceptibility curve for the 
range of heating rates used. 

Figure 1 shows the real part of the susceptibility curves 
measured for three separate films with Fe thicknesses of 
1.5 ML, 1.75 ML and 2.0 ML (all ±0.1 ML), and heat- 
ing rates varying from 0.03 K/s to 0.70 K/s. (Note: the 
data for 1.5 ML Fe is the same data that is analysed in 
reflU). The absolute scale of the susceptibility is un- 
certain within a factor of about 2 because the magneto- 
optic Voigt parameter is not well known for these films. 
However, the relative scale for all traces on all plots is 
consistently calibrated to the absolute optical rotation. 
The general shape of the curves is consistent with many 
previous studies. At high temperature the susceptibility 
decreases exponentially as the domain density increases. 
At low temperature the susceptibility decreases exponen- 
tially as the domain walls become pinned in structural 
defects. The peak in the susceptibility occurs at an in- 
termediate temperature due to the interrelation of these 
two processes, and depends strongly on the film thick- 
ness. 

The dependence of the susceptibility on the heating 
rate that is evident in fig. 1 is summarized more clearly 
in fig. 2 using two quantities: the temperature at the 
peak of the susceptibility, T p k, and the full width at half 
maximum of the peak. The systematic variation of these 
quantities with heating rate is the subject of section IV. 
The remainder of this section is concerned with compar- 
ing the variation of the susceptibility with film thick- 
ness to previous findings. Quantitative fits to the data, 
using eq.© with parameters defined in eq.Q and ©, 
are shown in fig. 3. Part a) demonstrates the quality 
of the fits and confirms that the essential points of the 
model are valid. Values of the parameters E p i n and n 
are given in parts b) and c), respectively. As has been 
shown previously^ the average pinning energy, which is 
a structural property of each film, is independent of the 
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FIG. 1: Magnetic ac susceptibility of x ML Fe/2 ML 
Ni/W(110) films as a function of temperature for different 
constant heating rates. All the measurements were performed 
after cooling the sample from 360 K at a rate of -0.10 K/s. 
An ac field of 2.0 Oe at a frequency of 210 Hz was applied 
normal to the film. Parts a), b) and c) show sequential mea- 
surements from a single film of thicknesses 1.5 ML, 1.75 ML 
and 2.0 ML of Fe, respectively. 



heating rate. The pinning energy decreases rapidly with 
film thickness. A previous study^ has shown that the 
pinning is due to changes in the perpendicular anisotropy 
caused by changes in the thickness at monolayer steps in 
the film, and that E p i n ~ d~ 3 < 2 for measurements of a 
single film grown sequentially to a number of thicknesses. 
The current results are qualitatively consistent with this 
finding, but quantitative comparisons are not possible 
among a series of independently grown films with unre- 
lated microstructure. The parameter n depends on both 
film thickness and the heating rate. The variation with 
thickness that displaces the curves one from another in 
fig. 3c is in qualitative agreement with eq.Q. Because 
the surface anisotropy varies as 1 /d, k can be ultimately 
derived from eq.([I|) as/tw ~ d ~ 3/2 - 

The thickness dependence of E P i n and k leads to the 
systematic displacement of the susceptibility curves with 
thickness observed in fig. 2. A gross measure of T p k is the 
pinning temperature T pini and eq.® indicates that this 
scales with E pin . Thus the curves in fig. 2a are displaced 
to lower temperature for thicker films. Similarly, since 
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FIG. 2: a) The temperature at the peak of the susceptibil- 
ity curve, T p k, as a function of the heating rate, for films 
of different Fe thickness. The arrows mark approximately the 
heating rate at which T p k has a minimum. The fitted lines are 
discussed in section IV. b) The full-width at half-maximum 
(FWHM) of the susceptibility curves as a function of heating 
rate. Lines simply connect the points. 



FIG. 3: Fits to the data using the phenomenological param- 
eters from eq.©. a) Representative fits for measurements at 
a heating rate of 0.30 K/s for different film thicknesses, b) 
The pinning energy for local Barkhausen steps of microscopic 
segments of a domain wall, c) The paramter re, which reflects 
the exponential increase in domain density with temperature. 
Lines simply connect the points. 



the two exponential factors (E p i n and k) that cut off 
the susceptibility above and below its peak value vary 
inversely with a power of the thickness, the peak width 
in fig. 2b increases as the film thickness increases. 



IV. ANALYSIS AND DISCUSSION 

As the heating rate increases, T p k first decreases 
sharply, and then reverses and gradually increases. The 
value of R at which this reversal occurs (marked by the 
arrows in fig. 2a depends systematically on the film thick- 
ness. This suggests that there are two competing dynam- 
ical processes - one that dominates at small R and the 
other that dominates at large R. It is important to note 
that the changes in the peak amplitude with R in fig. 1 
are very modest even as the peak temperature changes 
substantially. 

The rapid decrease in T p j. and in the peak width at 
small R are correlated to an increase in k in fig. 3. This 
behaviour was analysed in ref. [l5l where it was related to 



the slow relaxation of topological defects in the domain 
pattern after the sample is quenched from a high temper- 
ature, delocalized phase, as in section IIB. The present 
data now confirms this observation at other thicknesses. 
The second, less dramatic, process at large R is more ob- 
vious in the new data at thicknesses of 1.75 and 2.00 ML. 
According to the analysis in section IIC, the increase in 
T p k at large values of R is qualitatively consistent with 
the dynamics of domain density changes and growth on 
an intermediate time scale. This agreement suggests a 
quantitative analysis using a combination of eq.® and 

T pk (R) = T - A cxp(-< e/ f (R)) + BR. (18) 

Figure 4 presents the least-squares residuals of the fit 
of the experimental peak temperatures to this expres- 
sion by varying To, A, Ed and B for a range of values of 
the time constant rod- A consistent optimum value of 
log 10 Tod = 0.00 ± 0.05 is found. Table I gives the best fit 
parameters for each film thickness. Error estimates on 
these parameters, as well as on rod, are derived by hold- 
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FIG. 4: The sum of squared residuals for the independent fits 
of eg. (|18[) to the peak temperatures as a function of heating 
rate in fig. 2a. The residuals for the film with 2.0 ML Fe are 
multiplied by 10 to place them on the same plot. The residuals 
for this film are smaller because there are data points at fewer 
heating rates. 



Fe ML 


E d (K) 


To (K) 


A(K) 


B fit (s) 


B calc. (s) 


1.50 


1560±25 


271.8±0.6 


26.0±0.9 


2.7±1.8 


3.95±0.5 


1.75 


1390±13 


247.7±0.5 


33.8±0.8 


11.0±1.5 


8.91±0.5 


2.00 


1270±12 


234. 1±). 5 


38.1±0.7 


21.0±1.5 


16.1±0.5 



TABLE I: Parameters for fitting the data in fig. 2a to eq. (|18[) 
appear to the left of the double line. In all cases log 10 (rod) = 
0.00 ± 0.05. The parameter B calc. to the right of the double 
line is determined in fig. 5. 



ing all the other parameters constant, and finding the 
range of variation that changes the squared residuals by 
1.0. The fitted curves are superimposed upon the data 
in fig. 2a. 

These results confirm the earlier analysis of the data 
for the 1.5 ML Fe films^ where the term BR due to to 
the dynamics of domain density changes was neglected. 
Since B is found to be small for this thickness, the val- 
ues of the other parameters are essentially unchanged 
from those found in the previous study. The new data 
for thicker films show that, upon quenching from high 



temperature, the relaxation parameter rod and Ed for 
topological defects in the domain pattern do not depend 
strongly on the film thickness. It continues to take many 
minutes to an hour to remove these defects and for the 
magnetic susceptibility to take the equilibrium curve rep- 
resenting the ordered stripe phase. The fact that Ed does 
not scale with E p i n further substantiates the conclusion 
that the removal of the topological pattern defects is not 
limited by microscopic pinning mechanisms, but rather 
by the low probability that weak dipole interactions will 
drive the co-ordinated domain wall fluctuations required 
to make the mesoscopic changes required to remove the 
defect from the pattern. For all thicknesses, the defect 
relaxation is reflected as well in a decrease in k. As was 
previously noted^ the larger value of k in the presence 
of topological defects may represent their influence on the 
non-equilibrium free energy, creating an increase in the 
magnetic "stiffness" as measured by the susceptibility. 

Because the constant B is not small for the 1.75 and 
2.00 ML Fe films, the data cannot be reasonably fit with- 
out this term. This proves that even as the topological 
defects are relaxing, a second, faster relaxation mecha- 
nism is present in these films. The sign and order of 
magnitude of the second effect is consistent with the es- 
timates from eq. (jTTJ) for relaxation of the domain density, 
and indicate that the domain density and topological de- 
fect density are not strongly coupled, but relax rather 
independently. 

The interpretation of the constant B can be tested 
quantitatively. Fig. 5 gives an example calculation for 
the film with 1.75 ML Fe. The parameters A, k,tq p and 
Epi n are those determined in fig. 3a by the fit to eq.© 
for the data with i?=0.30. In fig. 5a the result for x eg (T) 
derived from A and k is given by the dashed line. The 
same parameters are then used to evaluate n eq (T) and 
derive n(T) as solutions to eq. (fT3"|) for different heating 
rates. The initial condition for each is the saturated, 
"frozen" domain density that results from solving eq. (|13[) 
as a function of cooling from the equilibrium state at 
350 K, at a rate of i?=-0.10 K/s. Once n(T) is known, 
Xeff(T) is easily derived from eg . (|T2"j) . The solid curves 
in fig. 5a give \eff(T) for heating rates of 0.20, 0.40, 
and 0.70 K/s (moving from left to right in the figure), 
for the choice a = 10 5 5 . It is clear that Xeff(T) peaks 
at lower temperature when R is smaller, since the system 
has more time to relax toward the equilibrium state. Fig. 
5b presents x(T) by multiplying the results of part a) 
by the dynamical pre- factor defined in eq.© and (jTTJ) . 
It can be seen that the calculated susceptibility shifts 
gradually to higher temperature with higher R, but that 
the amplitude of the peak changes very little. 

The position T p k and peak amplitude can be taken 
from a series of curves such as those in fig. 5b. The 
linear dependence of these quantities as a function of R 
can be extracted to give one point on each of fig. 5c and 
5d. Repeating the calculation for many values of a pro- 
duces the entire plots. The three traces on these plots 
correspond to parameters fit to the data in fig. 3a for 
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Simulation for 1.75 ML Fe/ 2.0 ML Ni/W(1 10) 
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imental data. These constraints determine the value of 
a for these films. The corresponding calculated values of 
B taken from fig. 5c are appended to Table I, under the 
heading "B calc". It can be seen that this single value 
of the single adjustable parameter a gives good quantita- 
tive agreement with "B fit" (within limits of uncertainty) 
determined from the independent analysis of the plots of 
T pk in fig. 2. 

These results are strong support for the identification 
of the dynamics of the domain density as the cause of 
the gradual rise in T p k at large heating rates. These find- 
ings validate the model presented in appendix A, where 
the domain density changes by domain growth driven by 
dipole interactions, but ultimately hindered by the same 
actviated domain wall pinning as the local movement of 
domain walls in response to an applied field. The value 
of a ~ 10 5 5 derived from the analysis demonstrates that 
the domain density relaxes at an intermediate time scale, 
that is 10 5 ' 5 times slower than the time scale for local pin- 
ning of a segment of a domain wall, but much faster than 
that required for the relaxation of topological defects in 
the domain pattern. 



V. CONCLUSIONS 



FIG. 5: Calculation of the changes in the susceptibility peak 
due to domain growth dynamics for a film with 1.75 ML Fe. a) 
Parameters fit in fig. 2 for the film with a heating rate of 0.30 
K/s are used to calculate Xe<j directly (dashed line), and Xeff 
through the integration of eq. (|13|) (solid lines). The heating 
rates are (left to right) 0.20, 0.40 and 0.70 K/s. The single 
adjustable parameter a = 10 5 ' 5 . b) The results in part a) are 
multiplied by the dynamic prefactor in eq.© to simulate the 
measured susceptibility, c) Calculations of the susceptibility 
as a function of R, such as those in part b), are used to find 
the linear dependence of T p k on R as as function of a. d) 
Calculations of the susceptibility as a function of R, such as 
those in part b), are used to find the relative change in the 
peak amplitude as a function of a. 



thicknesses of 1.5, 1.75 and 2.0 ML Fe. As was previ- 
ously discussed, for small values of a, the characteristic 
temperature for the freezing of the domain density, T„, 
is much less than that for the pinning of the domain wall 
segments, T p i n , and since any changes in domain density 
dynamics with R cannot be observed, B « 0. For a j 
10 5 , T p i n << T n , and T p k depends upon R through a 
relatively constant, non-zero linear coefficient B. In the 
intermediate range of a, T p i n ss T n and the two dynamic 
mechanisms interact in a complex manner that gives a 
strong sensitivity of B to the value of a. 

The two vertical lines on fig. 5c and 5d highlight the re- 
gion near a = 10 5 5 . In this region, B is positive and rela- 
tively insensitive to a, and the peak height depends only 
weakly on the heating rate, in agreement with the exper- 



Magnetic susceptibility measurements of Fe/2 ML 
Ni/W(110) films show the effects of three separate mag- 
netic relaxation mechanisms that span a wide range of 
fundamental time scales. Two of these are readily ap- 
parent and have been studied previously. Field-driven 
oscillatory movement of existing domain walls has a fun- 
damental time scale of 10 -9 s and is pinned by structural 
defects in the film. This dynamical response cuts off the 
susceptibility below a pinning temperature T p i n and gives 
the ac susceptibility its characteristic shape. Topological 
domain pattern defects that persist after quenching from 
a delocalized phase are removed on a fundamental time 
scale of 10° s, and cause large (25 to 40 K) shifts in the 
entire susceptibility curve when it is measured at a slow 
heating rate (R < 0.10 K/s). This shift is to lower tem- 
perature as R is increased. Through a systematic study 
of these effects, it has been possible in the present ar- 
ticle to isolate and study a third relaxation mechanism 
at an intermediate time scale. This is the dipole-driven 
changes in domain density that create a smaller shift of 
the ac susceptibility curve to higher temperature as R is 
increased (R > 0.20 K/s). The increase in characteris- 
tic times of the mechanisms reflects the increase of the 
size of feature upon which each mechanism works - from 
the magnetic response of individual Barkhausen steps, to 
the correlated fluctuations required to drive the growth 
of domains, to the even more complex rearrangements 
needed to resolve topological defects involving multiple 
domains. 

A simple activated relaxation model of domain den- 
sity dynamics provides a quantitative description of 22 
susceptibility traces measured at different heating rates 
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on films with three different thicknesses, using a single 
value of a single adjustable parameter. This parameter 
is a, the proportionality between the relaxation time for 
local domain wall pinning and domain density changes, 
T n (T) = aTp(T). The functional form of this relation 
gives strong support for three features of the model. Most 
fundamentally, it confirms that the changes in domain 
density in the experiment are driven by the departure of 
n(T) from n eq (T) in the absence of a magnetic field, that 
is by long range dipole forces. Second, that the same 
structural defects that pin local domain wall movement 
limit the growth of domains. Finally, that the nucleation 
of domains is not a limiting factor in changing the do- 
main density in these experiments, and can be neglected. 
This could be because the nucleation energy is less than 
the pinning energy, or due to the magnetic history of the 
films. Since all films are grown at high temperature, they 
start out with a large density of domains. Upon cooling, 
the domain density decreases by shrinking the domains, 
but may leave behind small nuclei which are not elimi- 
nated by the weak dipole forces.— Upon heating, growth 
can occur at the existing nuclei. 

The quantitative value a ~ 10 5 ' 5 confirms that the do- 
main density relaxes on an intermediate time scale. The 
separation of the three relaxation mechanisms by orders 
of magnitude permits them to be studied independently 
by appropriately designed experiments. In particular, 
changes in the domain density are found to occur inde- 
pendently from the resolution of topological defects, both 
because of the very different time scales, and the unre- 
lated activation energies. This statement does not mean 
that structural defects do not affect the local evolution of 
domain walls in the topological defects, but rather that 
they are not the limiting factor in their dynamics. These 
findings are in essential agreement with microscopy stud- 
ies. The very long relaxation times for topological defects 
permits them to be imaged, and spatial Fourier trans- 
forms of images show a loss of orientational order. The 
much quicker changes in stripe density, however, cannot 
be followed, except through the sudden disapperance of 
individual stripe domains between image line-scans i&i 3 - 

The present experiments have many similarities to sim- 
ulations by Bromley et al.£2- These authors also identify 
three dynamical regimes with different time scales for the 
relaxation of the magnetic state of Ising spins on a square 
lattice. There are also important differences, since the 
simulations study relaxation from a magnetically satu- 
rated state at constant temperature, rather than quench- 
ing and heating at (almost) zero field. Nonetheless, the 
initial, fast relaxation of the magnetization in the sim- 
ulations, driven by Zeeman energies once the field is re- 
moved, is clearly analogous to the field-driven domain 
wall oscillation in the current experiments. Similarly, the 
subsequent growth and connection of domain segments to 
form stripes on an intermediate time scale is very remi- 
nescent of the model used in the present article for the 
change the domain density. In the simulations, this pro- 
cess is driven both by dipole interactions and a residual 



non-zero magnetization. Finally, Bromley et alM> find 
that topological defects form boundaries separating re- 
gions of different stripe orientation, and require a very 
long time to relax. This is qualitatively similar to the 
relaxation of topological defects observed indirectly here. 
A comparative analysis should not be pressed too far, 
since the simulations are for very narrow Ising domain 
walls and do not include a mechanism for domain wall 
pinning. 
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Appendix A 

In this appendix, a relaxation equation for changes in 
the domain density is developed. Assume that the long 
domains characteristic of the stripe phase are aligned 
along the y-axis, and that the measurement samples an 
area of dimension L x x L y . If the i th domain segment has 
length £j, then the average domain density is the total 
domain perimeter divided by the sample area: 

u=-i-^2£ i . (Al) 

y 1 

The rate of change of the domain density is then 

— = — - = N 2vi. (A2) 

dt LxLy dt L x Ly 

Vi is the velocity of the growing (or contracting) tip of 
the domain segment. The factor of 2 enters because the 
length of a domain segment changes at both ends. 
Due to pinning, the activated average velocity is 29 

v = v cxp( e—^ — £ ). (A3) 

The fundamental speed |i>o| is given by the product of 
the Barkhausen step size s and the "attempt" frequency 
i>q = 1/tq p . The structural defects that pin the domain 
walls and retard the growth of the domain segments are 
the same ones that pin the domain walls and retard the 
change in the width of the domains when an ac field is 
applied. Therefore, E p i n and tq p are the same quantities 
in eq. dU) and (|A3I) . Given that these pinning energies are 
of order 10 3 K, and that the dipole energy Edi po i e driving 
domain formation is of order 10° K, the latter can be ne- 
glected in eq. (|A3I) . To the extent that the distribution of 
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pinning energies can be represented by an average energy 
E pin , then the distribution of domain wall velocities Mj 
can be represented by the average v. 

The sum L 1 L J2i 2 represents the density of active 
growth fronts. In a relaxation approximation, the density 
of active growth fronts is proportional to the total density 
of growth fronts, p, times the deviation of the number of 
domain walls from its equilibrium value, (N — N eq ). The 
total density of growth fronts includes the ends of existing 
domain segments as well as latent domains nucleated at 
defects that may grow if the conditions are favourable. 
This model does not include an explicit nucleation en- 
ergy, but rather assumes that, upon cooling, the weak 
dipole energy differences can reduce the domains centred 
on the nucleation centres to a size such that the domain 
perimeter is negligible compared to that of the remaining 
domain segments, but cannot completely remove them. 
This assumption is supported by microsopy studies of 
domain formation^ 

Putting these factors together, 



Equation[T3lis substituted in for §Zk , recalling the relation 



T n (T) = otT p (T). This gives 



2Epi n Rot 



r*(T pk ) 



u 1 +uj 2 r p l (T pk ) n(T pk ) 



(B2) 



As was discussed in section IIC, T pk will depend on the 
relaxation of the domain density only when T n > T pin . 
In this case uj 2 Tp in (T p k) « 1 and can be neglected in 
the denominator of the dynamical factor. A further con- 
sequence is that near T pk , n is substantially less than n eq , 
as can be seen in fig. 5a and 5b. Using these approxima- 
tions gives as result that is not valid in the limit of very 
small R: 



2E pin Ra 3 3/™ s 

-UJ T p (l pk ) 



■i(T p k) 



kT pk 



n(T pk ) 



(B3) 



dn 
~dt 



2(ispL x 



exp(- 



-E n 



kT 



-)(n - n eq ). 



(A4) 



The proportionality constant, f3, between the number of 
active growth sites and p(N — N eq ), is difficult to derive. 
It is related to how strongly the driving dipole energy 
depends upon the number of domain walls. To the ex- 
tent that dipole energies are "small" and depend only 
logarithmically on the number of domains, f3 is expected 
to be "small" and produce a long relaxation time. The 
factor \/(2pL x ) = £, the average separation of growth 
sites in the y direction. Finally, the differential equation 
describing the relaxation of the domain density is 



dn(T) 



-1 



where 



dt T n (T) 

otTp, with a = £/(s/3). 



(n(T) - n eq (T)), 



(A5) 



Appendix B 

To find how the temperature of the peak of the sus- 
ceptibility, Tpk, depends upon the heating rate, it is first 
necessary to find an expression for T p k- Starting from 
eq. (TTTT) . the derivative ^ is set to zero to yield the im- 
plicit relation 



2E, 



pin 



rl(T P k) 



1 



dn(T) 



kTp k 



+ u 2 T 2 (T pk ) n( T pk ) dT 



(Bl) 



To find the implicit derivative of the left hand 

side of this equation is straight forward. For the right 
hand side, 



d ( n eq (T pk ) 
dR \ n(Tp k ) 



1 



x I re 



n eq (T pk ) 
n(T pk ) 

dn(T) 



i(T pk ) dT 



T=T„. 



(ff ). (B4) 



where the ansatz in eq.((3]) has been used. 

Substituting in this equation from eq. dBip and (|B3I) . 
and setting the derivatives of the left and right hand side 
of eq. (|B3l) equal, yields after rearrangement 



i« ( ^k)(A 

y dR ' { T pk 



pk 



+ ^(l-^ 2 r 2 p (T pk )) + K )R 

2^2 



(B5) 



Using once again the condition u 2 Tp in (T p k) « 1 gives 
the final estimate 



dT 



pk 



dR 



(B6) 
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